Gamma-Z box contributions to parity violating elastic e-p scattering 
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Parity-violating (PV) elastic electron-proton scattering measures Q-weak for the proton, Q^, . To extract Q^, 
from data, all radiative corrections must be well-known. Recently, disagreement on the yZ-box contribution to 
2w has prompted the need for further analysis of this term. Here, we support one choice of a debated factor, 
go beyond the previously assumed equality of electromagnetic and yZ structure functions, and find an analytic 
result for one of the yZ-box integrals. Our numerical evaluation of the yZ-box is in agreement within errors with 
previous reports, albeit somewhat larger in central value, and is within the uncertainty requirements of current 
experiments. 
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Parity-violating (PV) elastic electron-proton scattering 
measures Q-weak for the proton, Q^, which is the Z-current 
Dirac form factor for the proton. This form factor has been 
measured at a number of energy scales, and part of the moti- 
vation for the PV experiment is to check its evolution against 
predictions of the standard model. Deviations from the pre- 
dictions could be a signal of new physics, that is, of currently 
unknown terms in the Lagrangian. Another motive is to check 
the consistency of Q!^ measured by NuTeV at low momen- 
tum transfer yj] versus measurements using other processes, 
although this motive may today be weaker due to the recogni- 
tion yfl of charge symmetry violating effects upon the NuTeV 
and other experiments with targets using large nuclei. 

Knowing all corrections is important to accurately obtain- 
ing Q'w from the e-p parity violating asymmetry. Correc- 
tions from, among other sources, yy, WW, and ZZ boxes 
have been well considered. A surprise came when Gorchtein 
and Horowitz J3| evaluated the inelastic corrections to the yZ 
box at zero overall momentum transfer (inelastic meaning the 
hadron state between the y and Z connections is not a proton), 
using a dispersive method that connected the box evaluation 
to the inelastic structure functions. Their result was unexpect- 
edly large and of of uncertain robustness. Sibirtsev et al. J3l 
subsequently reevaluated the yZ box, finding even larger cor- 
rections but importantly asserting that the uncertainty in the 
corrections was safely below the projected uncertainty in the 
experimental result. This is in line with the conclusions in fl. 

Because of the importance of the consequences and of the 
differences in the two results, including an overall factor "2" 
in one of the main formulas, we present another visit to this 
subject. We corroborate the factor "2" as given by Sibirtsev et 
al. and give numerical evaluations with at least partly different 
input that leads to results slightly larger but compatible within 
uncertainty limits compared to Sibirtsev et al. We also show a 
useful technical advance not mentioned in J^t] or 0], that one 
of the triple integrals required to obtain the answer can be ana- 
lytically done, leading to an easier numerical evaluation of the 
final result. Our results are still at zero momentum transfer; 
a partonic calculation of the yZ box valid at high momentum 
transfer can be found in ffj ■ 

Calculation of yZ box diagrams. The quantity is the 
Z-boson current Dirac form factor of the proton, evaluated at 



zero momentum transfer. One can measure it from the parity 
violating asymmetry in elastic electron-proton scattering, 



A PV = 



(1) 



where Or,l are cross sections for electron helicities A = ±1/2 
and unpolarized protons. To lowest order, Apv comes from 
interference between single y and single Z exchanges, and 
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where t is the overall momentum transfer, negative for space- 
like momentum transfers, and Qw = 1 — 4 sin 2 9w{0). With 
corrections, one has, following |7||, 
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The WW and ZZ box diagrams give Dww ar, d Dzz, and are 
well calculated perturbatively. The yZ box diagrams, Fig.Q] 
involve low momentum scales where perturbation theory is 
not reliable for the hadronic part of the diagram. Gorchtein 
and Horowitz J^, calculating only contributions from the 
inelastic intermediate states (elastic contributions have been 
considered in JI- 11]), showed how to dispersively relate the 
yZ box at t = to hadronic structure functions. With some 
approximations, they obtained a result that was larger than ex- 
pected. Sibirtsev et al. [Hf] improved the calculation, obtaining 
in fact a somewhat larger result but with tighter uncertainty 
limits. 

Demands on the uncertainty limits are set by current and 
planned experiments. The Q-weak experiment at JLab aims to 
measure Q 1 ^ to about 4% combined statistical and systematic 





FIG. 1: The y-Z box diagrams. 
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error at an incoming electron energy of 1 . 165 GeV and MAMI 
is discussing an experiment with 180 MeV incoming electron 
energy measuring to perhaps a part in a thousand. 

Theoretically, the quantity Dyz is obtained from the par- 
ity violating part of Jftyz, Fig.Q] by comparison to the corre- 
sponding term in single Z-exchange, 



□ r z = 



fyZ;A = l/2 " 



iyZ;X=-l/2 p ,LO 
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where W% = (M + m n ) 2 , m n is the mass of the pion, and 
Qmax = ( s —M 2 )(s — W 2 )/s. The Fj Z structure function has 
been dropped because the ratio gy/g e A ~ 0. The remaining 
term is labeled by a superscript "V" to indicated its association 
with the vector part of the Z-boson current. This expression 
agrees with the result reported in 0]. 

The real part, Re Oyz, is given by the dispersion relation, 



-^Z;A=l/2 - y ^Z;X=-\/2 

(Reference 13] presents results using 8yz = Oyz/Qw L ° ) 

For the exchange of a Z-boson between an electron with 
momentum k and proton with momentum p, the denominator 
of Oyz is 



Ren; z (£) = ^ 
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where V K = (W 2 —M 2 )/2M. 

Rewriting Eq.dTob in the current context as 



*Z;A = l/2- 
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A = 



(2ME') 2 + 2ME'(M 2 -W 2 - Q 2 ) -M 2 Q 2 
Q 2 {-M 2 + W 2 + Q 2 ) 



(13) 



where g e A = -|. 

The calculation of the numerator of Re Dy z requires the ap- 
plication of the optical theorem. The imaginary portion of the 
amplitude for photon, Z-boson exchange is 



we notice that the E' dependent terms can be separated from 
those dependent on Q 2 and W 2 . The E' integrands can be 
evaluated analytically if the order of integration is switched. 
Sliding the energy integration to the inside changes all three 
integration bounds. The W 2 integral is now evaluated from 
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q2 + ^ W 2_ M 2 + g2)2 +4M 2g2)l/2] tQ ^ After eva l uating 

the E' integrals, Re D^ z becomes 
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Here, k\ is the intermediate 4-momenta of the electron, Q 2 
—q 2 , and g e v = — j + 2sin 2 0w'- One obtains 



Emin+E J Q 2 (W 2 -M 2 + Q 2 ) 

Evaluation o/ReD^. Experimental data do not exist for 
7 ^2- In the scaling region, high Q 2 and high W, there are 



yz 

separated parton distributions 012111311 and one gets F-[ 9 using 
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similar to the purely electromagnetic F'i, where 2e q g q —> e q . 
However, one expects and can verify that the bulk of the sup- 
port for the yZ box comes from the resonance region and from 
lower Q 2 . In order to proceed, earlier work accepted lUHl the 
equality Fj Z = F^l, which can be shown to be approximately 
true in certain regions and certain limits. We will investigate 
the equality and improve upon it. 

Our numerical evaluation of ReD^ z uses the Christy - 
Bosted fits lfl4ll in the resonance region (W < 2.5 GeV), the 
Capella et al. fits in the high-energy low-Q 2 region (W > 2.5 
GeV and Q 2 < 5 GeV 2 ), both of these with some modifica- 
tion, and used the CTEQ parton distributions CT10.00 ll 211 in 
the scaling region (W > 2.5 GeV and Q 2 > 5 GeV 2 ). 
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For resonance ph otop roduction and electroproduction, the 
parton model (e.g., IU5I0 shows how each amplitude depends 
on the quark charges. It is useful to note that later analysis 
indicated that two-quar k op erators play a small role in photo- 
production amplitudes 111 611 . The charges then can be changed 
to the Z-boson vector coupling parameters g% to compare res- 
onance contributions in F^ and F^. For any isospin-3/2 res- 
onance, the result is just a multiplication by ( 1 + Q^ L ), since 
only the Al = 1 currents contribute ||4j] . 

Other resonances are more complicated. For example, for 
the D\j( 1520), the A3/2 amplitude also scales like (1 )> 
but the A 1 12 amplitude has two contributions, one sharing a 
matrix element with A 3 1 2 and one which will be multiplied by 



(1/3 + Q^ L0 ). There is extra phenomenological information, 
that the A3 ^ dominates in photoproduction and that there is 
a rapid transition to the high Q 2 dominance of the A 1 n ex- 



pected from hadron helicity conservation |17 



1811 . This give 

enough information to modify the Christy-Bosted contri- 
bution for the yZ in a Q 2 dependent fashion. As a remark, 
the average Q 2 within the integrals for incoming energies in 
the JLab range is only about 0.4 GeV 2 . Similar considera- 
tions apply to the 7^15(1690), although now the multiplication 
factors are Q^ L0 for the A 3/2 and (2/3 + Q'^ LO ) for the other 
amplitude, so the reduction from the purely electromagnetic 
case is quite noticeable. The modification of the resonant part 
of the Christy-Bosted fit is thus straightforwardly done, and 
gives a resonant contribution to Fj\ about 9% smaller than 
to FjT. The Christy-Bosted fits come within 3% of nearly all 
the data points, and the points themselves have comparable 
(mostly systematic) error. We allow some margin, assigning a 
10% uncertainty in this part of the calculation. 

An additional note is that the amplitude for electromagnetic 
excitation of a proton to a state with quark spin-3/2 is propor- 
tional to (e u + 2e c i) (this is the Moorhouse selection rule 11911 ). 
which is not zero when turned into its Z-current analog. How- 
ever, this excitation seems small also for a neutron target, so 
we do not consider it further. 

The resonance region fit includes a smooth background 
non-resonant part, which one can think of as scattering off 
collections of quarks with scant final state interactions. In 
a full SUf(3) limit, where all light quarks are equally likely 
and which may be pertinent in a high-energy xQ 2 / (2Mv) 
limit, one has f/f/^/I = 1 + Q^ L °- In a valence quark limit 

with SU(6) wave functions, one gets (2/3 + Q^ W ) for the 
same ratio. The latter is better at high-x and the former is bet- 
ter at low-x and we take the mean, and use the extremes to set 
our uncertainty estimate. One can examine the fJ^/fJ^ ra- 
tio in the scaling region, and the result along one boundary of 
the CTEQ region is shown in Fig. [2] The value at the low W 
end is in agreement with out high-x expectation for the back- 
ground in the resonance region. The rest of this CTEQ-based 
curve is at lower-x and the F^ and FjJ, structure functions are 
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nearly within 5% of equality for much of the range. This also 
marks what we may expect at the upper end of the Capella et 
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FIG. 2: The ratio FJ Z /FT* vs. W obtained from the CTEQ parton 
distribution functions at fixed Q 2 =5 GeV. 



al. region, with expectation of closer equality as Q 2 further 
decreases. We estimate the modification of the Capella et al. 
fit for the present case by multiplying it by a W dependent 
function which is the mean of unity and this boundary curve, 
and take the extremes to estimate the uncertainty. 

Our numerical results for ReD^ as a function incoming 
electron energy are shown in Fig. [5| We also show the sep- 
arate results from the resonance region and from above the 
resonance region, and show uncertainty limits for the total. 
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FIG. 3: Plot of Re vs. incoming electron lab energy. The dashed 
red line gives the resonance region contribution using the Christy 
and Bosted lfl4ll structure function fit; the dotted blue cur ve g ives 
the non-resonance region contribution using Capella et al. luOll and 
CTEQ et al. (TJ. The solid black curve gives the total, with an error 
band indicated. 
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However, the goal is now a part in a thousand measurement 
of Q^,, so the uncertainty requirement is about 0.00006. This 
remains a challenge, and will require further accurate fits in 
the resonance region, further thinking about the validity of the 
F^ = F\\ approximation, and considering the Re Q^ z term. 

We conclude by restating that the yZ box contribution is 
known well enough for current experiments, but that more ac- 
curate determinations will be needed within several years. 
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FIG. 4: Plot of F 2 {Q 2 l W) vs. W using (H at Q 2 = 1.525 GeV 2 . 



The data is from JLab 121112 



For the JLab Q-weak experiment, Ei a b = 1.165 GeV, and 

ReD^ z (1.165 GeV) = 0.0057 ±0.0009. (16) 
This agrees within uncertainty limits with the Sibirtsev et al. 
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For information, at the JLab energy, the high Q 2 above the 
resonance region contribution from CTEQ gives 0.00019 of 
the total; this is about 0.00004 lower would be gotten by sim- 
ply extending Capella et al. to the high Q 2 region. Also, 
the longitudinal part of the structure functions contribute only 
about 0.0007 to the above result, roughly evenly split among 
resonances in the resonance region, non-resonant background 
with W < 2.5 GeV, and contributions where W > 2.5 GeV. 

Though we agree with the Sibirtsev et al. result, part of the 
agreement is due to the reduction in our result from analyzing 
the F 2 = FP relation. Had we used the equality everywhere 
but the scaling region, our result would have been 0.00065 
higher. One difference between us is that in the resonance re- 
gion, we used the Christy-Bosted fit lfl4ll . which represents the 
data to 3% or better over almost the entire applicable range. 
By way of examples, Christy and Bosted give plots of cross 
section vs. W at a number of incoming energies and angles. 
Sibirtsev et al. used their own dedicated resonance region fits, 
and also fit the data well, as seen in their plots of F2 vs. W at 
several fixed Q 2 's fl]. To facilitate direct comparison, Fig. [4] 
here shows the Christy-Bosted F2 vs. Wata typical Q 2 . 



Discussion. Using the Q 1 = value 
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0.23867(16) fl2j then Q% w = 0.04532(64) and with the cor- 
rections listed in Eq. d3}, 2w somewhat exceeds 0.07. A 4% 
measurement of this number requires an absolute accuracy of 
about 0.0028, so the ReD^ z correction needs to be known 
more accurately than, say, 30%. This we believe is the case. 

For the future, discussions of a PV experiment at lower en- 
ergy 112411 are partly prompted by the smaller expected yZ box 
correction. For 180 MeV we obtain 
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ReDy Z (180 MeV) = 0.00125 ± 0.00018 . (17) 



